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Abstract. We present high spatial resolution 
(~0'/6) near-infrared broad-band JHK images and 
Br7 2.1661 /im and H2 1-0 S(l) 2.122 /im emission line 
images of the nuclear regions in the interacting starburst 
galaxies NGC 520, NGC 1614 and NGC 7714. The 
near-infrared emission line and radio morphologies are in 
general agreement, although there are differences in de- 
tails. In NGC 1614, we detect a nuclear double structure 
in BT7, in agreement with the radio double structure. We 
derive average extinctions of Ak = 0.41 and Ak = 0.18 
toward the nuclear regions of NGC 1614 and NGC 7714, 
respectively. For NGC 520, the extinction is much higher, 
Ak — 1.2 - 1.6. The observed H2/Br7 ratios indicate 
that the main excitation mechanism of the molecular gas 
is fluorescence by intense UV radiation from clusters of 
hot young stars, while shock excitation can be ruled out. 

The starburst regions in all galaxies exhibit small 
Bi - 7 equivalent widths. Assuming a constant star forma- 
tion model, even with a lowered upper mass cutoff of M M 
= 30 Mq, results in rather old ages (10 - 40 Myr), in 
disagreement with the clumpy near-infrared morpholo- 
gies. We prefer a model of an instantaneous burst of 
star formation with M„ = 100 M , occurring ^6-7 Myr 
ago, in agreement with previous determinations and with 
the detection of W-R features in NGC 1614 and NGC 
7714. Finally, we note a possible systematic difference in 
the amount of hot molecular gas between starburst and 
Seyfert galaxies. 
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1. Introduction 

Interaction and merging play an important role in the evo- 
lution of galaxies. Interacting galaxies usually have strong 
infrared (IR) emission (e.g. Rieke et al. 1980; Telesco, Wol- 
stencroft & Done 1988), and show a wide range of nuclear 
activity, including AGN, nuclear starbursts (SB), ultralu- 
minous IR galaxies and post-SB activity. These phenom- 
ena probably all arise from the ability of interactions to 
transport gas into the galactic nuclei (e.g. Barnes & Hern- 
quist 1991) to fuel both AGN and SBs. In SB galaxies, the 
radiation output is dominated by star formation (SF), and 
SBs typically have SF rates of 5-50 M Q yr _1 within a re- 
gion of 0.1-1 kpc extent, much larger than that found in 
the Galactic center (e.g. Mezger, Duschl & Zylka 1996) or 
in normal galaxies (e.g. Keel 1983), and they thus provide 
an extreme environment in which to study the SF process. 

Multiwavelength studies of interacting SB galaxies can 
provide insights into the connection between interactions, 
massive star clusters, and nuclear activity. The relatively 
unobscured near-IR (NIR) emission provides a better han- 
dle than optical emission to quantify the properties of SF 
in galaxies. In this paper we present high resolution NIR 
broad-band JHK images and Br7 and H2 1-0 S(l) emis- 
sion line images of the circumnuclear regions of three in- 
teracting/merging SB galaxies, NGC 520, NGC 1614 and 
NGC 7714. Because of their proximity, brightness and rea- 
sonably large angular size, they provide excellent targets 
for these studies. Br7 originates from H II regions sur- 
rounding hot young OB star clusters, while H2 arises from 
hot molecular gas and traces the material available for SF. 

In the remaining part of Section 1, we give brief intro- 
duction to the galaxies. In Section 2 we discuss the ob- 
servations and data reduction. In Section 3 we discuss the 
morphology of the circumnuclear regions, determine the 
extinctions to the SF complexes, constrain their SF prop- 
erties, stellar populations and SF history by comparison 
among the NIR tracers and with multiwavelength emis- 
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sion, and discuss the gas masses of the galaxies. Conclu- 
sions are drawn in Section 4. 

1.1. NGC 520 

NGC 520 (Arp 157, UGC 966, VV 231) at v sys = 2059 
km s- 1 (D = 27.3 Mpc, 1" = 130 pc, for H = 75 Mpc" 1 
km s _1 and qo = 0.5) is an IR-luminous peculiar pair 
of galaxies with a highly disturbed flattened morphology 
(e.g. Stanford & Balcells 1990; Hibbard & van Gorkom 
1996). The inclination of NGC 520 is -66° (Rownd & 
Young 1999). NGC 520 is at an intermediate merger stage 
(Hibbard & van Gorkom 1996), with the two nuclei em- 
bedded within a common stellar envelope. The primary 
nucleus (hereafter PN) is optically hidden behind a dust 
lane at PA = 95°, while the second, optically bright nu- 
cleus (hereafter NWN) is situated -40" (5.3 kpc) to NW. 
A bright optical tail stretches — 2.'8 (22 kpc) to SE, bend- 
ing sharply to E and connecting to a broad stellar plume. 

The PN has -4.3 x 10 9 M of molecular gas in a -7" 
x 3" (920 x 400 pc) scale rotating molecular gas disk at 
PA = 95° (Yun & Hibbard 2000). While the CO emission 
is in good agreement with the 1.4 GHz radio emission (5" 
x 2" = 660 x 260 pc, PA = 93°; Condon ct al. 1990), both 
avoid the Ha emission, which is dominated by plumes of 
ionized gas emerging up to a projected distance of 20" (2.6 
kpc) from the nucleus along PA = 25° (Hibbard & van 
Gorkom 1996). These plumes probably represent a SB- 
driven bipolar outflow of ionized gas out of the nuclear SB 
region (Norman et al. 1996). While Ha emission regions 
exist within the PN region, the nuclear SB is completely 
obscured optically by dust within the nuclear gas disk (but 
becomes visible in Br7 ; Section 3.1.1). 

NGC 520 is probably the result of an encounter —3 
x 10 8 yr ago between a gas-rich and a gas-poor galaxy 
(Stanford & Barcells 1991). Both the S tail and the plume 
arose from the gas-poor progenitor and are currently dis- 
turbing the underlying gas distribution from the gas-rich 
disk (Hibbard & van Gorkom 1996). The optical spec- 
trum of the NWN is dominated by A stars, indicating that 
it is now in a post-SB phase (Stanford & Balcells 1990; 
Bcrnlohr 1993b), in agreement with the non-detection of 
Br7 emission (Section 3.1.1). Enhanced SF is currently 
detected only in the PN (e.g. Stanford 1991). 

1.2. NGC 1614 

NGC 1614 (Mrk 617, Arp 186, II Zw 15) is an IR-luminous 
SB(s)c pec type SB galaxy at v sys = 4723 km s _1 (D = 
62.4 Mpc, 1" = 300 pc). It has been extensively studied at 
optical (e.g. Dc Robertis & Shaw 1988), NIR (e.g. Aitkcn, 
Roche & Phillips 1981; Forbes et al. 1992; Puxley & Brand 
1999) and radio (e.g. Condon et al. 1982) wavelengths. The 
inclination of NGC 1614 is -30° (Rownd & Young 1999). 

The extended emission around NGC 1614 is highly 
asymmetric, suggesting recent tidal interaction with at 



least one other galaxy and that the galaxies are cur- 
rently merging. The chaotic optical structure includes two 
roughly symmetrical inner spiral arms, a linear tail to the 
SW extending almost 1' (18 kpc) from the nucleus and 
a large curving arc to the E/SE that reaches —33" (9.8 
kpc) from the nucleus (e.g. Neff et al. 1990). Only a single 
peak is seen in the NIR continuum, suggesting that the two 
galaxies have already merged. However, two peaks (possi- 
bly double nuclei) are seen in the Br7 emission (Section 
3.1.2) and in radio (Condon et al. 1982). 

Armus, Hcckman & Miley (1989) found a weak broad 
feature at 4660 A due to W-R stars in the nuclear spec- 
trum, which was later confirmed by Vacca & Conti (1992). 
Neff ct al. (1990) made a multiwavelength study of NGC 
1614 and concluded that its extreme IR luminosity results 
from vigorous nuclear SF induced by the interaction and 
merger of at least two galaxies. The molecular hydrogen 
mass of NGC 1614 is M(H 2 ) = 1.1 x 10 10 M (Sanders, 
Scovillc & Soifer 1991). There is an unresolved central CO 
source in NGC 1614 which emits -30 % of the total CO 
flux (Scoville et al. 1989), implying a very high molecular 
gas concentration to the nucleus. 

1.3. NGC 7714 

NGC 7714 (Arp 284, VV 51, Mrk 538) is a nearby (v sys = 
2808 km s _1 , D = 37.2 Mpc, 1" = 180 pc) peculiar SBb 
type SB galaxy. It is very IR-luminous and has strong 
SB activity both in the nucleus and in circumnuclear re- 
gions (e.g. Gonzalez-Delgado et al. 1995, hereafter GD95; 
Garcia- Vargas et al. 1997). NGC 7714 and NGC 7715 (at 
-2' (21 kpc) to the E; v sys = 2770 km s"\ D = 36.7 
Mpc) form the Arp 284 system of interacting spirals, as 
is evident from the HI and optical bridge (Smith, Struck 
& Pogge 1997; Papaderos & Fricke 1998) between them. 
The inclination of NGC 7714 is -42° (Chapelon, Contini 
& Davoust 1999). A striking feature of NGC 7714 is an 
asymmetric circumnuclear ring of diameter —20" x 40" 
(3.6 x 7.2 kpc), with most of the emission located to the 
E of the nucleus. This ring is largely devoid of Ha (GD95) 
and HI emission (Smith et al. 1997), and has NIR colours 
similar to the old disk population (Bushouse & Werner 
1990). This structure is thus probably related to the dy- 
namical perturbation of the disk rather than to previous 
SF. A weak He II 4686 A W-R feature has been detected 
in NGC 7714 (GD95). The molecular gas mass of NGC 
7714 is M(H 2 ) = 2.1 x 10 9 M Q (Sanders et al. 1991). 

While NGC 7714 shows strong SB activity, NGC 7715 
has no optical emission lines and is in a post-SB phase 
(Bcrnlohr 1993a). There are three main circumnuclear HII 
regions in NGC 7714, labelled A, B and C by GD95, and 
located in the bulge/disk of the galaxy at 5" (900 pc) E, 
12" (2.1 kpc) NW and 14" (2.5 kpc) SW from the nucleus, 
respectively. Region A is detected by us as Br7 region 4. In 
the deconvolved Ha image of GD95, the nuclear SB region 
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breaks furthermore into the nucleus and an extranuclear 
component, which we identify as Br7 region 2. 

2. Observations and data reduction 

The observations of the Br7 2.1661 /jm and H2 1-0 S(l) 
2.121 /zm emission lines and the JHK bands were car- 
ried out in September 1998 with the 3.8 m United King- 
dom Infrared Telescope (UKIRT) on Mauna Kea, Hawaii, 
under FWHM 0"6-0"7 seeing. We used the 256x256 px 
IRCAM3 camera, with pixel size 0'.'281 px -1 and field of 
view ~ 70" x 70". For the emission line observations, we 
used cooled (T — 77 K) narrow-band filters and a Fabry- 
Perot (F-P) etalon with spectral resolution ^400 km s _1 
and equivalent width (EW) 0.0038 /im. For the emission 
lines, sequences of 1 minute jittered observations in a five 
point grid were made at the on-line wavelength and at the 
nearby blue and red continuum. For the JH K -bands, 30 
sec integrations were made both centered on the galaxy 
and on the adjacent sky. The total integration times were 
40, 10 and 30 minutes in Br7, 23, 10 and 30 minutes in 
H 2 and 2.5 - 7.5 minutes in the JH if-bands, for NGC 
520, NGC 1614 and NGC 7714, respectively. 

The line images were linearized, dark-subtracted, flat- 
fielded and sky-subtracted using IRAFQ These images 
were aligned to within a small fraction of a pixel using 
field stars or the nucleus as reference, and merged into one 
on-line image and two continuum images. The continuum 
images were scaled, combined, and subtracted from the 
on-line images, and the final line images were flux cali- 
brated against spectral type A standard stars. 

The sensitivity of the F-P depends on both the line- 
of-sight (l.o.s.) velocity of the galaxy and the shift in the 
transmitted wavelength of the F-P over the field of view. 
These effects were corrected for by dividing the measured 
fluxes by an inverse Airy function (e.g. Bland-Hawthorn 
1995). For NGC 520 and NGC 7714, we used the l.o.s. 
Ha velocity fields by Bernlohr (1993b; PA = 120° and PA 
= 145°) and GD95 (PA = 110° and PA = 216°), respec- 
tively. For NGC 1614, we are not aware of a published 
velocity field, and we used the Ha N-S rotation curve by 
De Robertis & Shaw (1988) assuming axisymmetry. In the 
nuclei of all galaxies, we assumed the l.o.s. rotational ve- 
locity to be km s . The emission line images have not 
been corrected for by the inverse Airy function, because 
noise would then depend on the position in the image. 
However, the observed emission line fluxes have been cor- 
rected for. To enhance S/N and to detect faint structures, 
all images were slightly smoothed to <1" resolution. 

Photometry in all bands was performed at the loca- 
tion of the detected Br7 emission regions. The aperture 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Fig. 1. The if-band image of NGC 520 with logarithmic 
intervals. The inset shows the central 6" x 6" of NGC 520 
PN. In this and all subsequent figures, north is up and east 
to the left. 

used for each region was selected to include as much emis- 
sion as possible, while avoiding overlap with neighbouring 
regions. The smallest distance between the nearest emis- 
sion regions in these galaxies is ~0"8, which is larger than 
the seeing during the observations (~0"6). The observed 
fluxes were corrected for Galactic extinction and redshift 
(K-correction). We estimate a photometric accuracy in 
the JHK magnitudes ~0.03 mag, in the JHK colours 
~0.05 mag, and in the emission line fluxes ~10 %. 

3. Results and discussion 

3.1. Morphology 
3.1.1. NGC 520 

The K-h&nd image of NGC 520 is shown in Fig. |] (see 
also e.g. Stanford & Balcells 1990). In the J-band (not 
shown), absorption of the stellar light distribution along 
the PN dust lane is noticeable (as in the optical) . However, 
a stellar bulge, bisected by the dust lane, begins to become 
evident in the PN. The overall appearance of the PN is 
that of an edge-on stellar disk. In the J-band, the NWN 
stellar bulge still dominates the emission from the NGC 
520 system. On the other hand, in the -RT-band, the dust 
lane no longer significantly obscures the underlying stellar 
light in the PN. The stellar bulge is seen clearly at the 
position of the center of the optical dust lane. The PN 
is elongated in roughly E-W direction, while the outer 
bulge region appears more spherical. The position of the 
-ftT-band PN coincides with the radio continuum and CO 
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Fig. 2. The Br7 emission of NGC 520 PN. The lowest con- 
tour is at 19 % level of the maximum and corresponds to 
3a. The other contours are at la intervals. The maximum 
surface brightness is 2.4 x 10~ 15 erg s _1 cm~ 2 arcsec -2 . 
The black dots correspond to the three Br7 emission re- 
gions studied in this work. 



peaks. The NWN is —4 times fainter than the PN in the 
i-C-band, indicating that it is the less massive of the two 
bulges (stellar masses of the old population are 9 x 10 9 
M and 2 x 10 9 M for the PN and NWN, respectively; 
Stanford & Balcells 1990). 

We have detected Br7 emission at higher than 3cr level 
in the central -5" x 3" (660 x 400 pc) region of NGC 
520 PN (Fig. ||). This emission describes a flattened mor- 
phology and breaks up into several components. We have 
selected three emission regions for further study, labelled 
E (east), C (center) and W (west). The total Br7 flux from 
NGC 520 PN, -1.2 x 10~ 14 erg s" 1 cnr 2 is in reasonable 
agreement with previous values, e.g. 1.9 x 10~ 14 erg s _1 
cm -2 in a 6" x 8" aperture (Stanford 1991). 

There is very good spatial correspondence between the 
brightest regions in the Br7 and the radio emission. At 

5 GHz, NGC 520 PN is an extended linear triple radio 
source with dimensions -6" x 1'.'8 (790 x 240 pc) at PA 
= 93° (Condon et al. 1982). At 15 GHz, the emission is 
resolved into five components at PA —90° (Carral, Turner 

6 Ho 1991). Because the non-thermal radio emission is 
believed to arise in supernova (SN) remnants and the ther- 
mal radio emission to be reradiated UV emission from OB 
stars, this correspondence is not surprising. 

The correspondence between the H2 (Fig. ||) and the 
Br7 emission in the PN is reasonably good. Both emission 
lines define an elongated structure in roughly E-W direc- 




Fig. 3. The H 2 emission of NGC 520 PN. The lowest con- 
tour is at 28 % level of the maximum and corresponds to 
3<7. The other contours are at la intervals. The maximum 
surface brightness is 0.93 x 10 -15 erg s _1 cm -2 arcsec -2 . 



tion. The H2 emission is more extended, with maximum 
dimensions of —8" x 4" (1.1 x 0.5 kpc). There is, however, 
no clear correlation between the Br7 and H2 peaks. The 
total H 2 flux from NGC 520 PN, -0.7 x 10~ 14 erg s" 1 
cm~ 2 is smaller than the flux 1.9 x 10 -14 erg s _1 cm~ 2 
in a 6" x 8" aperture (Stanford 1991). However, as the 
spectrum of Stanford (1991) covers only a few pixels, sug- 
gesting the possibility of contamination from alignment 
and/or bad pixels, this comparison has to be approached 
with caution. 

The H-K colour map of NGC 520 PN is shown in 
Fig. ^| The nuclear H-K colour is very red, corresponding 
to extinction of A^ —1.5. Note that the reddest peak is 
displaced from the if-band nucleus by — 0'.'3 (40 pc) to 
the E. The size and the PA of the red region is similar to 
that of the Br7 emission. Outward the extinction quickly 
becomes much less severe. There is no spatial correlation 
between the Ha (Bernlohr 1993b) and the Br7 emission 
in the PN, again indicating that the effect of extinction is 
severe in NGC 520 (see also Section 3.2.1). 

3.1.2. NGC 1614 

The X-band image of NGC 1614 (Fig. |) clearly shows the 
inner spiral arms, but is not deep enough to show the SW 
tail. However, the SF regions SW and W of the nucleus 
are detected, although they remained undetected in the 
NIR image of Neff et al. (1990). At very faint levels, we 
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Fig. 4. The H-K colour map of NGC 520 in greyscale and 
contours. The highest contour corresponds to H-K = 1.3 
and the other contours are at H-K = 0.15 intervals. 




Fig. 5. The if-band image of NGC 1614 with logarithmic 
intervals. The inset shows the central 4" x 4" region. 



also see evidence for the beginning of the optical curved 
arc to the E of the nucleus. 

We have detected Br7 emission at higher than 3cr level 
in the central ~3'/5 x 3'.'5 (1.0 x 1.0 kpc) region of NGC 
1614 (Fig. §). Intriguingly, this emission is resolved into 
two main components (SE and NW) straddling the nu- 
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Fig. 6. The Br7 emission of NGC 1614. The lowest con- 
tour is at 7.1 % level of the maximum and corresponds to 
3cr. The other contours are at 9.5, 14, 21, 31, 43, 57, 74 
and 86 % of the maximum. The maximum surface bright- 
ness is 7.5 x 10~ 15 erg s" 1 cm~ 2 arcscc~ 2 . The black dots 
correspond to the two Br7 emission regions studied in this 
work. 



cleus and separated by ~0'/7 (210 pc) in roughly SE-NW 
direction. The total Br7 flux from NGC 1614, -2.4 x 
10~ 14 erg s -1 cm -2 is somewhat smaller than previously 
reported values, e.g. 7.2 x 10~ 14 erg s _1 cm~ 2 in a 7" 
aperture (Ho, Beck & Turner 1990), 3.6 x 10" 14 erg s _1 



-2 
-14 



en) - in a 3" x 4" aperture (Neff et al. 1990) and 6.0 x 

erg s _1 cm~ 2 



10" i4 erg s" 1 citT 2 (Puxley & Brand 1999). This differ- 
ence suggests that a part of the more extended Br7 emis- 
sion may have remained undetected in our observation, 
due to the relatively short integration time. 

This is the first time that the double structure is seen 
in the NIR. Note that it is not visible in our (Fig. |J) or 
published if-band images (e.g. Neff et al. 1990; Forbes et 
al. 1992), due to poorer spatial resolution and heavy fore- 
ground extinction. However, their existence was suggested 
by Puxley & Brand (1999) who detected two compo- 
nents of roughly similar brightness in a velocity-resolved 
Br7 spectrum of the nucleus of NGC 1614. 

The 5 GHz radio emission of NGC 1614 (Neff et al. 
1990) is in good agreement with the Br7 emission, showing 
two similar brightness maxima NW and SE of the nucleus, 
separated by ~l'/2 (360 pc). The total extent of emission 
is -3" x 3" (890 x 890 pc) for both radio and Br 7 . The 
center of NGC 1614 is resolved into two nuclei of similar 
brightness also at mid-IR 12.5 /im wavelength (Keto et al. 
1992), oriented ~E-W and separated by ~0"8 (240 pc). 



6 



J.K. Kotilainen et al.: NIR line imaging of starburst galaxies 




2 0-2 
A RA(") 



Fig. 7. The H 2 emission of NGC 1614. The lowest contour 
is at 22 % level of the maximum and corresponds to 3cr. 
The other contours are at la intervals. The maximum 
surface brightness is 1.8 x 10~ 15 erg s _1 cm -2 arcsec -2 . 
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Fig. 8. The H-K colour map of NGC 1614. The highest 
contour corresponds to H-K ~ 0.64 and the other con- 
tours are at H-K — 0.08 intervals. 



An arm reaches N of the W component and then arches 
over to the E. This mid-IR structure resembles the double 
source observed in the 5 GHz radio emission and in Br7, 
although the orientation is slightly different. 

The double peaks in NGC 1614 are probably analogous 
to the double nuclei observed in many other IR-luminous 
galaxies. The separation of the nuclei in NGC 1614, ~0'/7 
(210 pc) is slightly smaller than e.g. in Arp 220 (330 pc; 
Graham et al. 1990), and considerably less than the sep- 
arations of 2-6 kpc observed in several ultraluminous IR 
galaxies (Carico et al. 1990). On the other hand, NGC 
1614 shows the structural characteristics of a collision with 
small impact parameter. The small separation of the nu- 
clei and the high central molecular gas density (Scoville 
et al. 1989) should be inducive for AGN-type activity. 
However, multiwavelength information (Neff et al. 1990) 
strongly suggests that the nuclear emission of NGC 1614 
is dominated by SB activity. 

The correspondence between the H 2 (Fig. [?]) and 
the Bi - 7 emission in NGC 1614 is reasonably good. The 
H 2 emission is slightly more extended than Bi'7, with 
maximum dimensions of —4" x 4" (1.2 x 1.2 kpc). The 
H2 emission is resolved into several peaks, but there is 
no clear correlation with the Br7 peaks. The H-K colour 
map of NGC 1614 is shown in Fig. [| The central region is 
much redder than the rest of the galaxy. The H-K map 
resolves this region of high extinction into a circunmuclear 
ring-like structure, with the Br7 maxima situated at the 
inner edges of this ring. 



3.1.3. NGC 7714 

The if -band image of NGC 7714 is shown in Fig. |[ In ad- 
dition to the strong nuclear emission and the inner spiral 
arms, there is some evidence for the detection of the ~20" 
x 40" (3.6 x 7.2 kpc) scale stellar ring to the E of the 
nucleus. We have detected Br7 emission at higher than 
3<7 level in the central ~6" x 7" (1.1 x 1.3 kpc) region 
of NGC 7714 (Fig. |Io| ). This emission is not symmetric, 
but consists of the nuclear emission and a chain of Br7 re- 
gions (2, 3 and 4) toward SE of the nucleus. The strongest 
Br7 region (1) is not coincident with the nucleus, but is 
situated ~0'/8 (140 pc) to SW. Its photometry was per- 
formed in a rectangular aperture to avoid the extended 
emission N of the nucleus and the nuclear emission. In- 
terestingly, there is also indication of a fainter peak ~0'.'4 
(70 pc) to NE of the nucleus, with spiral-like emission con- 
necting the two peaks. If real, this would add NGC 7714 
to the increasing number of active galaxies with nuclear 
SF rings. The total Br 7 flux of NGC 7714, -2.7 x 10~ 14 
erg s _1 cm -2 is slightly smaller than previously reported 
fluxes, e.g. 4.7 x 10~ 14 in a 7" aperture (Ho et al. 1990). 

The 5 GHz radio emission of NGC 7714 (Condon et 
al. 1982) shows two main components at PA —45° sepa- 
rated by ~l'/2 (210 pc), with the SW component being 
stronger. This morphology is in good agreement with the 
Br7 emission. There is a rather good correspondence be- 
tween the H 2 (Fig. 11) and the Br7 emission in NGC 7714. 
The maximum dimension of the detected H 2 emission is 
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Fig. 9. The if-band image of NGC 7714 with logarith- 
mic intervals. The inset shows the central 4" x 4" region. 
There is a foreground star toward SW of the nucleus. 
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Fig. 11. The H 2 emission of NGC 7714. The lowest con- 
tour is at 31 % level of the maximum and corresponds to 
3<7. The other contours are at la intervals. The maximum 
surface brightness is 5.9 x 10 -16 erg s _1 
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Fig. 10. The Br7 emission of NGC 7714. The lowest con- 
tour is at 4.1 % level of the maximum and corresponds to 
3cr. The other contours are at 5.4, 8.2, 12, 16, 22, 29, 46, 
57, 69, 83 and 98 % of the maximum. The maximum sur- 
face brightness is 7.0 x 10 -15 erg s _1 cm -2 arcsec -2 . The 
black dots correspond to the four Br7 emission regions 
studied in this work. 



~6" x 8" (1.1 x 1.4 kpc), and it is resolved into several 
peaks, in reasonable agreement with Br7. 

The nuclear SB extends to -16" x 12" (2.9 x 2.1 kpc) 
at PA— 110° in the Ha emission (GD95). This extended 
nuclear SB is surrounded by a circumnuclcar ring of small 
HII regions at —8" (1.4 kpc) radius and two giant HII 
regions at 12" (2.1 kpc) and 22" (3.9 kpc) from the nu- 
cleus. There is a good correlation between the Ha regions 
of GD95 and the Br7 emission, assuming that the nu- 
clear Ha emission is coincident with our region 1. Region 
4 also has a counterpart in the Ha image to within O'.'l. We 
have not detected Br7 emission further out from the nu- 
cleus, as opposed to the Ha imaging of GD95, suggesting 
a low value of extinction outside the nuclear region. The 
H-K colour map of NGC 7714 is shown in Fig. g The 
effects of extinction are the least severe of all the galaxies 
studied here. The reddest regions in the H-K map are in 
good agreement with the morphology of the Br7 emission 
(Fig.©- 

3.2. Extinction 

We have estimated the extinction towards the ionized 
sources in NGC 1614 and NGC 7714 from published re- 
combination line fluxes ( Aitken et al. 1981; Taniguchi et 
al. 1988; Armus et al. 1989; Puxley & Brand 1994; Storchi- 
Bergmann et al. 1995; Calzetti 1997). We have assumed 
standard interstellar dust properties, case B recombina- 
tion, and the interstellar extinction law Ay oc A~ 1,85 (Lan- 
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Fig. 12. The H-K colour map of NGC 7714. The high- 
est contour corresponds to H-K = 0.35 and the other 
contours are at H-K = 0.05 intervals. 



dini et al. 1984), and calculated the average integrated ex- 
tinction applicable to the central region of both galaxies. 
For NGC 520 PN, this method can not be used since it is 
optically obscured, and the extinction was assumed to be 
that calculated from the continuum colours (see below). 

Calzetti, Kinney & Storchi-Bergmann (1996) found 
that the geometry of the dust obscuring the ionized gas 
in a sample of 13 SBs can be well described by either 
homogeneous or clumpy foreground distribution and that 
the reddening is compatible with a Galactic-type extinc- 
tion curve. Thus, most of the dust is located around but 
outside the SB region. We are unable to justify either a 
homogenous or a clumpy foreground dust screen based on 
our data. For simplicity and for easy comparison with pre- 
vious results, we have assumed the case of a homogenous 
foreground dust screen, but for caveats, see Kotilainen et 
al. (2000). Note that the exact value of extinction does 
not affect the EWs of the emission lines, if the differential 
extinction between the lines and the continuum is small 
(e.g. Calzetti 1997). 

The extinction towards the continuum sources was es- 
timated by comparing the observed JHK colours of the 
emission regions to those of normal unobscured spiral 
galaxies (Glass & Moorwood 1985). This method, how- 
ever, may be biased toward low extinction regions, as the 
most heavily reddened regions may not be detectable in 
the JHK images. Also, the intrinsic colours of the galax- 
ies may be bluer than in normal galaxies, leading to an 
underestimation of the continuum extinction. 



3.2.1. NGC 520 

We used the continuum method for NGC 520, since 
Ha emission is completely obscured in the PN (Hibbard & 
van Gorkom 1996). We derive for the emission regions in 
NGC 520 PN very large extinctions of A K = 1.2-1.6, with 
the central component C having the largest extinction. 
Our extinction determination is in rather good agreement 
with previous results, Ay = 7-14, from [S III] /Ha ratio 
(Young, Kleinmann & Allen 1988), from Br7/Ha (Stan- 
ford 1991) and from NIR colours (Bushouse & Werner 
1990; Stanford & Balcells 1990). 

3.2.2. NGC 1614 

The central region of NGC 1614 suffers heavy dust ex- 
tinction, although less severe than in NGC 520 PN. We 
derive extinction of Ak = 0.41 from recombination line 
fluxes for the central region of NGC 1614. Very similar 
values are derived from the NIR colours, Ak = 0.41 and 
0.38 for the SE and NW regions, respectively. These val- 
ues are in good agreement with previous results, Ay = 3 
- 5, from hydrogen recombination line ratios (Puxley & 
Brand 1994; Shier, Rieke & Rieke 1996), from Pa/3/[Fe II] 
ratio (Puxley et al. 1994) and from NIR colours (Neff et 
al. 1990; Oliva et al. 1995; Shier et al. 1996). 

Puxley & Brand (1994) determined from NIR spec- 
troscopy that the dust distribution in NGC 1614 is best 
described either by clumpy foreground dust, implying Ay 
~4.7, or by homogeneous foreground dust mixed with in- 
ternal dust, implying Ay ^15. The fact that the levels 
of extinction deduced from the line ratios and the NIR 
colours are very similar indicates that a foreground screen 
of equal strength for stars and ionized gas is a good ap- 
proximation for NGC 1614 (see also e.g. Shier et al. 1996). 

3.2.3. NGC 7714 

The central region of NGC 7714 suffers the least amount 
of extinction of the three galaxies studied, with Ak = 
0.18 derived from the recombination line fluxes. Again, 
very similar values are derived from the NIR colours for 
the emission regions, 0.13 < Ak < 0.21 with, interestingly, 
the strong near-nuclear peak showing the smallest amount 
of extinction. Thus, the circumnuclear SBs may be more 
dusty and/or more evolved than the nuclear burst (see 
also Gonzalez-Delgado et al. 1999). The derived extinction 
values are in good agreement with previous results, Ay = 
1-2, from hydrogen recombination line ratios (Puxley & 
Brand 1994 ) and from NIR colours (Oliva et al. 1995). 

3.3. Star forming properties 

In Table |l] we give for the Br7 emission regions in all galax- 
ies the displacement from the if-band nucleus, the aper- 
ture diameter, the calibration coefficient from the Airy 
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Region 


ARA 


ADEC 


Apert 


corr a 


Br 7 


H 2 


J-H 


H-K 








tt 




10~ 15 ergs 


— 1 —2 

s cm 


mag 


mag 


NGC 520 E 


1.6 


0.1 


1.4 


1.11 


2.14 


0.965 


1.79 


1.19 


NGC 520 C 


0.1 


0.0 


1.4 


1.06 


2.58 


0.738 


1.80 


1.27 


NGC 520 W 


-1.5 


0.3 


1.4 


1.04 


1.98 


0.545 


1.72 


1.01 


NGC 520 ALL 


0.0 


0.0 


•5.6 




11.9 


7.29 


1.41 


1.00 


NGC 1614 SE 


0.2 


-0.3 


0.8 


1.05 


4.15 


0.816 


0.73 


0.48 


NGC 1614 NW 


-0.3 


0.3 


0.6 


1.11 


1.88 


0.264 


0.72 


0.47 


NGC 1614 ALL 


0.0 


0.0 


3.9 




23.8 


9.22 


0.74 


0.49 


NGC 7714 1 


-0.4 


-0.7 


2.8 x 1.7 


1.04 


11.9 


1.47 


0.63 


0.30 


NGC 7714 2 


2.6 


-0.6 


1.7 


1.01 


1.16 


0.443 


0.65 


0.34 


NGC 7714 3 


3.4 


-1.7 


1.1 


1.01 


0.346 


0.211 


0.67 


0.36 


NGC 7714 4 


4.5 


-2.7 


1.7 


1.01 


0.897 


0.171 


0.64 


0.34 


NGC 7714 N 


0.0 


0.0 


3.4 




17.5 


2.79 


0.63 


0.31 


NGC 7714 ALL 


0.0 


0.0 


11.2 




27.3 


8.35 


0.64 


0.28 



a Correction factor for the velocity field. 

Table 2. Dereddened Br 7 emission regions in the galaxies 



Region 


A a 


AO 
A K 


Bry 


H 2 


£?7"y 


J-H 


H-K 




mag 


mag 


10~ 15 ergs 


— 1 — 9 

s cm 




mag 


mag 


NGC 520 E 




1.50 


8.86 


4.23 


0.477 


0.13 


0.22 


NGC 520 C 




1.62 


12.0 


3.64 


0.303 


0.00 


0.22 


NGC 520 W 




1.21 


6.22 


1.79 


0.288 


0.38 


0.22 


NGC 520 ALL 




1.21 


37.4 


24.0 


0.641 


0.07 


0.22 


NGC 1614 SE 


0.41 


0.41 


6.12 


1.22 


0.200 


0.28 


0.22 


NGC 1614 NW 


0.41 


0.38 


2.80 


0.395 


0.141 


0.27 


0.20 


NGC 1614 ALL 


0.41 


0.41 


35.2 


13.8 


0.392 


0.28 


0.22 


NGC 7714 1 


0.18 


0.13 


14.1 


1.76 


0.125 


0.43 


0.19 


NGC 7714 2 


0.18 


0.19 


1.38 


0.529 


0.384 


0.45 


0.22 


NGC 7714 3 


0.18 


0.21 


0.410 


0.252 


0.615 


0.47 


0.24 


NGC 7714 4 


0.18 


0.18 


1.06 


0.204 


0.192 


0.44 


0.22 


NGC 7714 N 


0.18 


0.14 


20.8 


3.33 


0.160 


0.43 


0.20 


NGC 7714 ALL 


0.18 


0.10 


32.4 


9.96 


0.308 


0.44 


0.17 



a Determined from literature emission line fluxes, assuming A\ cc A -1,85 (Landini et al. 1984). 
Estimated from the H-K colour. 



function, and the observed Br 7 and H2 fluxes and JHK 
colours. In Table || we give the extinction derived from 
the recombination line fluxes and from the H-K colour, 
and the dereddened Br 7 and H2 fluxes, H2/Br 7 ratio and 
JHK colours for the Br 7 regions. The smallest values of 
FWHM of these regions, - 1'.'4 correspond to a size 

of the emitting region of 150 - 250 pc. Therefore, the emis- 
sion regions detected in the NIR actually are conglomer- 
ates of several OB associations and giant molecular clouds, 
probably similar to scaled-up versions of the 30 Dor H II 
region in the LMC (e.g. Walborn et al. 1999). 

The H2/Br 7 ratio can give important clues about the 
excitation mechanism(s) of the hot molecular gas (Pux- 
ley, Hawarden & Mountain 1990). The main mechanisms 
suggested are thermal excitation in hot gas by low veloc- 
ity shocks (e.g. Draine, Roberge & Dalgarno 1983) or by 



intense X-ray radiation (e.g. Maloney, Hollenbach & Tie- 
lens 1996), and fluorescent excitation by strong UV radi- 
ation (e.g. Black & van Dishoeck 1987). The dereddened 
H 2 /Br 7 ratios span a range 0.29 - 0.48 in NGC 520 PN, 
0.14 - 0.20 in NGC 1614 and 0.12 - 0.62 in NGC 7714. 
None of these ratios can be explained by shock excita- 
tion caused by cloud-cloud collisions or a SN-driven wind 
(H2/Br 7 > 1). Interestingly, only two of the Br 7 regions 
(NGC 520 E and NGC 7714 3) can readily be explained 
with fluorescent UV excitation by individual hot young 
stars (H2/Br 7 = 0.4 - 0.9), while for the large majority of 
them, the line ratios of < 0.4 are in agreement with flu- 
orescent excitation by intense UV radiation from a large 
compact cluster of hot stars (Puxley et al. 1990). 

For NGC 1614, our line ratios are in good agreement 
with those derived by Goldader et al. (1997; 0.214 in a 3" 
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x 12" aperture) and Moorwood & Oliva (1988, 0.203 in a 
6" aperture). Our ratios for NGC 7714 agree well with that 
by Taniguchi et al. (1988; 0.172 in a 3"5 x 7" aperture) 
and the upper limit by Moorwood & Oliva (1988; < 0.333 
in a 6" aperture). 

The SF properties were interpreted using the com- 
prehensive stellar population synthesis model for galax- 
ies with active SF by Leitherer et al. (1999; hereafter 
L99). It implements the latest stellar evolutionary tracks 
(Charbonnel et al. 1999) and atmosphere models (Leje- 
une, Buser & Cuisinier 1997), but does not treat self- 
consistently chemical evolution, binary evolution, mass 
loss and mixing processes, and late phases of stellar evolu- 
tion (L99). Despite these shortcomings, it may be consid- 
ered the most up to date of the suite of models available 
for studying SF galaxies (see e.g. Leitherer et al. 1996 and 
references therein). In particular, although the exact age 
of the SB is rather model-dependent, the general conclu- 
sion of young vs. old SB remains model-independent (see 
discussion in Reunanen et al. 2000). 

The L99 model predicts the evolution of NIR, optical 
and UV spectral features as a function of the burst age, 
metallicity, and the initial mass function (IMF) with lower 
and upper mass cutoff and slope a, for the limiting cases 
of instantaneous burst of SF (ISF) and constant SF rate 
(CSFR). The key parameter is the EW of Br7, which is 
sensitive to the age of the SB since it measures the ra- 
tio of young hot stars (Br7) and evolved RSGs (if-band 
continuum) independently of extinction. The estimation 
of EW(Br7) was unfortunately complicated by the im- 
possibility to adequately subtract a de Vaucouleurs bulge 
model (depicting the old stellar population; Kotilainen et 
al. 2000) from the if -band image in these disturbed, pe- 
culiar galaxies. The Br7 EW was thus determined by di- 
viding the Br7 fluxes with the total -RT-band fluxes. The 
resulting EWs are, therefore, lower limits and the derived 
SB ages upper limits only. However, for strong compact 
SBs as in these galaxies, a large fraction of the continuum 
light probably arises from the SB population and there- 
fore the derived lower and upper limits are probably not 
far from the real values. 

The L99 models predict the number of ionizing pho- 
tons below 912 A, N(H°), which can also be estimated 
from the Br7 flux. N(H°) allows us to evaluate the mass 
of recently formed stars (in ISF) or the SFR via an as- 
sumed IMF (in CSFR). Since we do not have enough data 
for more detailed modelling, and to allow for an easy com- 
parison with previous results, we assume solar metallicity, 
a = 2.35, and consider two models in what follows: (1) 
ISF with M u = 100 M and (2) CSFR with M u = 30 M Q . 
The SF properties (N(H°), EW(Br 7 ), age, mass, SFR and 
usn) of the emission regions in the galaxies were deter- 
mined by comparing the observed quantities with the L99 
models, and are given in Table |. 

The lower limits for the Br7 EWs are small, 11-14 A, 
9-10 A, and 20 - 39 A, for NGC 520 PN, NGC 1614 and 



NGC 7714, respectively. Assuming the CSFR model with 
M u = 30 Mq, we derive relatively high upper limits for 
the burst ages, 17-26 Myr, 36-38 Myr, and 9-11 Myr. 
With the exception of NGC 7714, these ages get close to 
the lifetime of individual giant molecular clouds in our 
Galaxy, after which turbulence and heating from SNe dis- 
rupt them and inhibit further SF (-20-40 Myr; e.g. Blitz 
1991). These ages are also in disagreement with the very 
clumpy NIR morphologies of the galaxies. The correspond- 
ing SFR is 0.9 - 1.8 M yr" 1 , 2.1-4.8 M Q yr" 1 , and 0.1 
-3.9 M yr" 1 , mass is 18 - 47 x 10 6 M Q , 8 - 17 x 10 7 
M Q , and 1.1 - 44 x 10 6 Mq, and SN rate is 9-25 xl0~ 3 
yr" 1 , 0.037-0.087 yr" 1 and 0.4-19 xlO -3 yr" 1 per region 
for the three galaxies. 

Assuming the ISF model, we derive much shorter up- 
per limits for the ages, 6-7 Myr for NGC 520 PN and 
NGC 1614, and -6 Myr for NGC 7714. The correspond- 
ing masses of the emission regions are 28-49 xlO 6 M , 
7.5-16 xlO 7 M , and 2 - 100 xlO 7 M , SFR is 4.4 - 7.5 
M Q yr" 1 , 11 - 25 M yr" 1 , and 0.4 - 16 M Q yr" 1 , and the 
SN rate is V SN = 26-46 xlO" 3 yr" 1 , 69-150 xl0~ 3 yr" 1 , 
and 2.4-97 xl0~ 3 yr" 1 , for NGC 520 PN, NGC 1614 and 
NGC 7714, respectively. 

The total mass of the hot gas in the Br7 regions of the 
galaxies, assuming the ISF model, —1.2 x 10 s M , —2.4 
x 10 8 M and ~1.2 x 10 8 M , although large, is still only 
a small fraction of the total molecular gas mass, —4.3 x 
10 9 M (Yun & Hibbard 2000), -1.1 x 10 10 M (Sanders 
et al. 1991) and -2.1 x 10 9 M (Sanders et al. 1991), in 
NGC 520 PN, NGC 1614 and NGC 7714, respectively. 

Stanford (1991) derived from extinction-corrected 
Ha fluxes SFR = 0.7 M yr" 1 for NGC 520 PN. The 
SFR in the PN is much higher than the SFR for a normal 
isolated spiral galaxy (0.02 M yr -1 ; Keel 1983). On the 
other hand, the SFR derived from the bolometric luminos- 
ity is 7 M yr" 1 for NGC 520 PN (Stanford 1991). Thus, 
the SB in the PN is unable by an order of magnitude to 
produce the bolometric luminosity, possibly caused by the 
old stellar population accounting for some fraction of the 
luminosity, or because the SB has just ended. 

The EW(Br 7 ), the EW(CO) and the N(Lyc)/L(bol) 
ratio all indicate a burst age between 6 and 8 Myr for NGC 
1614 (Puxley & Brand 1999), implying that the most mas- 
sive stars have already disappeared. The CO spectroscopic 
index = 0.26 indicates that about half of the if-band light 
arises from the old population. For the case of ISF, Puxley 
& Brand (1999) derive the mass of formed stars to be 0.3 
- 25 x 10 9 M for different IMFs. The very low M/L ratio 
of NGC 1614 (0.003; Joseph & Wright 1985) also excludes 
a stellar population older than —16 Myr. The conclusion 
of a short duration (< 1 Myr), recent (< 6 Myr) burst of 
SF is strongly supported by the presence of W-R features 
from a population of massive luminous young stars (Vacca 
& Conti 1992). 

Garcia- Vargas et al. (1997) modeled the NGC 7714 
circumnuclear SB regions with a young burst of age 3.5 - 
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Table 3. Star formation properties of the galaxies. 



instantaneous star formation constant star formation rate 



Region 


N(H°) 


EW(min) 


age 


mass 


SFR a 


VSN 


age 


mass 6 


SFR 


VSJV 




10 52 s" 1 


A 


Myr 


10 6 M Q 


M Q yr" 1 


10" 3 yr -1 


Myr 


10 6 M 


M Q yr" 1 


10" 3 yr -1 


NGC 520 E 


6.03 


14.1 


6.46 


40.3 


6.23 


36.8 


16.6 


21.8 


1.31 


11.1 


NGC 520 C 


8.16 


11.4 


6.55 


49.0 


7.47 


46.0 


26.4 


46.8 


1.77 


24.8 


NGC 520 W 


4.23 


13.2 


6.47 


28.5 


4.40 


25.9 


19.4 


17.8 


0.92 


9.38 


NGC 520 ALL 


25.5 


9.20 


6.61 


142. 


21.4 


131. 


40.5 


162. 


4.01 


71.6 


NGC 1614 SE 


21.7 


9.72 


6.59 


165. 


25.0 


151. 


36.0 


171. 


4.75 


87.0 


NGC 1614 NW 


9.96 


9.45 


6.60 


75.4 


11.4 


69.5 


38.4 


82.6 


2.15 


37.5 


NGC 1614 ALL 


125. 


9.35 


6.61 


949. 


144. 


876. 


39.1 


1050. 


26.9 


471. 


NGC 7714 1 


17.8 


19.9 


6.36 


103. 


16.2 


97.5 


11.4 


44.5 


3.90 


18.5 


NGC 7714 2 


1.74 


22.0 


6.33 


9.53 


1.51 


9.01 


10.8 


4.11 


0.381 


1.61 


NGC 7714 3 


0.518 


28.0 


6.25 


2.47 


0.40 


2.37 


9.82 


1.12 


0.114 


0.39 


NGC 7714 4 


1.34 


39.1 


6.13 


5.42 


0.88 


4.86 


9.02 


2.68 


0.297 


0.82 


NGC 7714 N 


26.3 


15.3 


6.43 


144. 


22.5 


134. 


14.2 


68.6 


4.83 


33.3 


NGC 7714 ALL 


40.9 


12.6 


6.50 


238. 


36.7 


214. 


21.3 


160. 


7.50 


84.3 



a : Mass divided by age 
b : SFR multiplied by age 



5 Myr to explain the emission line spectrum and the de- 
tected W-R features (GD95). In the optical spectrum of 
their region A (our region 4) 5" SE of the nucleus, they de- 
tect Ca II triplet absorption 8600 A features, compatible 
with two populations, the ionizing population with age ~5 
Myr, and a relatively young RSG-rich population of ~10 
Myr, responsible for the Ca triplet and the Balmer ab- 
sorption. Further modeling of the UV - NIR spectrum of 
the nucleus of NGC 7714 (Gonzalez-Delgado et al. 1999) 
found a best fit for a —4.5 Myr burst, with upper mass 
cutoff > 40 M Q . Gonzalez-Delgado et al. (1999) derive a 
SN rate = 0.007 yr -1 both from the spectral modelling 
and from radio emission, consistent with the nondetection 
of SNe in NGC 7714 in the monitoring by Richmond, Fil- 
ippenko & Galisky (1998). 

The M/L ratio of NGC 7714 is 0.020 (Bernlohr 1993a), 
while normal spirals span a range of 0.3 < M/L < 2.8. The 
small M/L ratios found in all three galaxies indicate an 
IMF biased toward high mass stars, probably through an 
increased lower mass cutoff (see also e.g. Rieke et al. 1980, 
Wright et al. 1988). The low M/L ratios also indicate that 
the SB has already produced enough RSGs to dominate 
the NIR emission over the RGB stars from the old popula- 
tion, indicating that the mass of the SB population must 
be at least a few per cent of the total mass. 

3.4- Gas masses 

Using the average surface brightnesses of the nuclear H2 
emission of the galaxies in the largest aperture (see Ta- 
ble ||), assuming T,,^ = 2000 K, and following the method 
of Meaburn et al. (1998), we derive -1000 M , -3000 
M Q and —800 M Q for the mass of the excited nuclear 



H 2 in NGC 520 PN, NGC 1614 and NGC 7714, respec- 
tively. These values should be multiplied by an unknown 
but probably small factor for the linewidth, which may 
be broader than the width of the F-P passband (see Sec- 
tion 2). The resulting masses are similar to those derived 
for the SB NGC 7771 (1700 M ; Reunanen et al. 2000) 
and the SB/Seyfert NGC 3079 (1200 M Q ; Meaburn et 
al. 1998), but much larger than derived for the Seyferts 
NGC 1097 (120 M ; Kotilainen et al. 2000), NGC 6574 
(80 M Q ; Kotilainen et al. 2000) and NGC 3227 (400 M ; 
Fernandez et al. 1999). Interestingly, in this small sample 
of galaxies, the SBs have a much larger H2 mass than the 
Seyferts. This comparison should be applied to a much 
larger sample of SBs and Seyferts to verify any difference 
in the amount of hot molecular material, and to study its 
implications for the fuelling of nuclear activity. 



4. Conclusions 

We present high spatial resolution (— 0'.'6) near-infrared 
broad-band JHK images and Br7 2.1661 /Ltm and H2 1-0 
S(l) 2.122 /im emission line images of the nuclear regions 
in the interacting starburst galaxies NGC 520, NGC 1614 
and NGC 7714. The near-infrared emission line and radio 
morphologies are in general agreement, although there are 
differences in details. In NGC 1614, we detect a nuclear 
double structure in Bi'7, in agreement with the radio dou- 
ble structure. We derive average extinctions of Ak = 0.41 
and Ak = 0.18 toward the nuclear regions of NGC 1614 
and NGC 7714, respectively. For NGC 520, the extinction 
is much higher, Ak = 1.2 - 1.6. The observed Bi2/Br7 ra- 
tios indicate that the main excitation mechanism of the 
molecular gas is fluorescence by intense UV radiation from 
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clusters of hot young stars, while shock excitation can be 
ruled out. 

The starburst regions in all galaxies exhibit small 
Br7 equivalent widths. Assuming a constant star forma- 
tion model, even with a lowered upper mass cutoff of M u 
= 30 M Q , results in rather old ages (10 - 40 Myr), in dis- 
agreement with the clumpy near-infrared morphologies. 
We prefer a model of an instantaneous burst of star forma- 
tion with M u = 100 M Q occurring ~6-7 Myr ago, in agree- 
ment with previous determinations and with the detection 
of W-R features in NGC 1614 and NGC 7714. Finally, we 
note a possible systematic difference in the amount of hot 
molecular gas between starburst and Seyfert galaxies. 
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